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The implementation of a three-level Lambda System in artificial atoms would 
allow to perform advanced control tasks typical of quantum optics in the solid state 
realm, with photons in the /im/mni range. However hardware constraints put an ob¬ 
stacle since protection from decoherence is often conflicting with efficient coupling to 
external fields. We address the problem of performing conventional STImulated Raman 
Adiabatic Passage (STIRAP) in the presence of low-frequency noise. We propose two 
strategies to defeat decoherence, based on “optimal symmetry breaking” and dynamical 
decoupling. We suggest how to apply to the different implementations of superconduct¬ 
ing artificial atoms, stressing the key role of non-Markovianity. 
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1. INTRODUCTION 


In recent years several experiments have demonstrated multilevel coherence in 
superconducting artificial atoms, as the observation of the Autler-Townes (AT) | l]|2j 
effect, of electromagnetically induced transparency (EIT) [3j, besides evidences of 
three-state superpositions [4j and coherent population trapping (CPT) (5J. Further 
exploiting coherence in such systems would be important in principle and moreover 
allow important applications in solid-state quantum integrated coherent architectures. 
So far all the experiments in these systems (except the one of Ref. |’5|) have been 
performed driving by ac-fields in ladder configuration (see Fig. |T]i). In this work 
we address the design of a lambda configuration in three-level artificial atoms which 
would allow to implement tasks [[6]-|8j where two-photon absorption and emission 
are invoved at once. Despite of several theoretical proposal |[9]-jT3|. this goal is still 
experimentally unsettled, mainly because protection from low-frequency noise re¬ 
quires to enforce exact or approximate symmetries of the Hamiltonian, which on the 
other hand imply selection rules canceling the pump coupling |[8|[9||~i~2| (see Fig.jTJi). 

Our second goal of is to elucidate the central role of non-Markovian noise in 
producing three-level decoherence for the class of phenomena based on CPT. We 
focus on a protocol called STIRAP [14] 151, described in Sec]2j which involves sev¬ 
eral basic coherent effects and allows striking applications in integrated atom-cavity 
systems. Therefore its demonstration would be a benchmark for multilevel advanced 
control in artificial atoms. In Sec. [3] we introduce an effective model for noise and 
argue that dephasing in the “trapped subspace” span{|0), 11)} (see Fig.|TJi) plays the 
major role. We show that implementation of STIRAP in Lambda configuration is 
possible within present technology. In Sec. [4] we propose two strategies to defeat de¬ 
phasing, namely the search for optimal symmetry breaking conditions, and selective 
dynamical decoupling of noise sources achieved by operating on a specific external 
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Fig. 1 - (a) Three-level system driven with AC fields in A configuration (in the insets the Ladder and the 
Vee configurations), (b) Gaussian pulses in the counterintuitive sequence (here flgT = 20, r = 0.6 T). 
(c) Population histories Pi(t) = |(i|-t/>(f))| 2 for ideal STIRAP (<5 = 0) and for 5 P = — 0.2flo> k= 1. 

control. Both strategies leverage on the fact that dephasing in the solid state is due to 
broad band colored noise (BBCN), which is inherently non-Markovian. As a conse¬ 
quence BBCN impacts on dephasing in a way specific of correlations of the induced 
fluctuations of the device bandstructure. Finally in Sec j6] we conclude and discuss 
some further perspective. 


2. COHERENT POPULATION TRANSFER IN THREE-LEVEL ATOMS 


STIRAP is an advanced control technique of M > 2-level systems, allowing 
complete population transfer between two states |0) and jf), even in absence of a 
direct coupling, via one or more intermediate states which are never populated. In 
three-level systems the indirect linkage is provided by the typical configurations of 
two ac-fields shown in Fig. [TJi. The pump field at u p « \Eo — Eq\, triggers transitions 
10) FA 12) whereas the Stokes, u s ~ \E >2 — E\ |, triggers |1) fa |2) ones. The standard 
Hamiltonian in the rotating wave approximation (RWA) in a rotating frame referred 
to the “bare” basis {10), 11), 12)} is given by the matrix 


H = 


0 
0 


0 

m 


S p (t) 


(i) 


where the Rabi frequencies 17 *. (t) for k = p,s are related to the amplitudes of the 
pump and Stokes fields, 6k are the single-photon detunings and 6 = <) p — d iS is the two- 
photon detuning. We will mostly refer to the Lambda configurations where 6 p (t) := 
E '2 — Pa — Up and 6 s (t) := E 2 — E\ — u s . At two-photon resonance, 5 = 0, the 
Hamiltonian (jT|) has an instantaneous eigenvector with null eigenvalue, eo = 0, given 
by | D) = (H s |0) — Q p \ 1))/+ ftp- It is called the “dark state” since state |2) 
is not populated, despite of the transitions triggered by the fields. In ideal STIRAP 
(5 = 0) adiabatic pulses are shined in the counterintuitive sequence, i.e. the 
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Fig. 2 - Left panel: typical Landau-Zener (LZ) patterns as a function of time of instantaneous eigen¬ 
values for nonzero detunings. Three different patterns may occurr according to the value of 5 p /5 (see 
Sec|4j. Right panel: efficiency diagram of STIRAP vs detunings. The curves enclose regions corre¬ 
sponding to efficiency larger than 90%. The black line is obtained for k p = k s and shows that the 
protocol is much more sensitive to deviations of <5 rather than of 8 P . Gray lines are the efficiencies 
when k p = 2 k s and k s = 2k p (see Sec. [4j. Intercepts of such curves with lines S p /5 = a define the 

two-photon linewidths 5i as a function of all the parameters but <5. The dashed line shows the example 

2 

of the Cooper pair box for q g = 0.48 and J — 1. 


Stokes preceding the pump. We will make use of Gaussian pulses 

n p = K P n 0 e ~ Kt+T)/T]2 Ms = K S n 0 e -^~ T/T]2 


( 2 ) 


with r ~ T. Here Ho is a frequency scale and k/- ~ 1 are constants which will be 
taken equal to 1 when not otherwise specified. In this way the dark state /.)(()) 
(Fig. [ljf) performs the desired |0) —> |1) evolution, yielding complete population 
transfer, while |2) is never populated (Fig. |T|l). 

Adiabaticity in ideal STIRAP [ 14JT5J requires that Hq T > 10. Since it involves 


in a clever sequence several coherent phenomena [ 151 (AT effect, EIT and adiabatic 
passage), STIRAP is very efficient, faithful and stable apart for a crucial sensitivity 
to 5 (see Fig. [2]right panel). Indeed for O' / 0 no exact dark state exists providing the 
adiabatic connection 10) -A 11). Still population transfer may take place by non-ideal 
STIRAP , via non-adiabatic Landau-Zener (LZ) transitions between adiabatic states 
(see Lig. [2] left panel), a mechanism crucial for the applications in artificial atoms 
(SecQ. 

In artificial atoms the Hamiltonian reads [11 16] H = F/o(q) + [A v (t) cos(uj p t) + 
A s (t)cos(u s t)]V, where the device Ho depends on tunable parameters q. The field 
couples to the system operator V and the envelopes Ak(t),k =p,s are slowly vary¬ 
ing with respect to Rabi frequencies. Under suitable conditions, H can be truncated 
to three levels. Performing the RWA and transforming to a doubly rotating frame, we 
get the form (jTJ, where Q p (t) = A p (t)V 2 o and Q s (t) = A s (t)V-n- 
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3. EFFECTIVE MODEL FOR SOLID-STATE NOISE 


Physically, noise in solid-state devices has large low-frequency components 
with a 1 //“ spectrum, and high frequency component, either white or ohmic. As¬ 
suming for simplicity a single noise source inducing fluctuations of the parameter 
q g , we can describe this BBCN by the phenomenological Hamiltonian [ 161 H = 
Ho(q g +x(t)) +A{t)V + H env . Here x(t) is a classical stochastic process accounting 
for low-frequency noise, whereas H env describes an environment coupled to the sys¬ 
tem, responsible for Markovian quantum noise. The effect of low-frequency noise is 
obtained by averaging over the stochastic process the density matrix p* (t\q g + x(t)), 
accounting for fast noise in a background stochastc field 


Pit) = JT>x(t) P[x(t)\ p f (t\q g +x(t)) 


( 3 ) 


Leading effects are estimated by evaluating the integral in the “quasistatic” or static 
path approximation (SPA), i.e. by substituting x(t) with a random variable x with 
distribution p(x) and calculating pf ( t\q g + x) by a Markovian master equation. 

Notice that Hq, its eigenenergies Ej(q g ) and the matrix elements Pj g enter¬ 
ing Qij(q g ) depend on q g . A proper choice of q g may enforce symmetries of H, 
which protect the system against dephasing due to fluctuations of Ei{q g ), but at 
the same time suppress some Pij. Non-Markovian noise determines fluctuations 
of the entries of the Hamiltonian (jTJ, namely 5k(q g + x) = A E k (q g + x) — ay- and 
VLk(q g + x) = Ak{t)Vk(q g + x) , where AE k (q g + x) and V k (q g + x) are the relevant 
energy splittings and ’’dipole” matrix elements. This is a key issue for all our sub¬ 
sequent analysis about design and optimization of Lambda systems. For instance, it 
is clear that for a Lambda configuration at nominal resonance, i.e. if external fields 
are resonant at the nominal bias q g , fluctuations in the “trapped subspace” translate 
in stray < 5/0 which are the most detrimental for STIRAP. It is convenient to expand 
detunings and Rabi frequencies. For instance at nominal resonance and for small 
enough fluctuations, imposing Eq = 0 we have 

5(x) = Ai(q)x+^Hi(q)x 2 ; 5 p (x) = AE 2 = A 2 {q)x+l ) B 2 {q)x 2 (4) 

where A;(q) = dE % (q)/ dq g and B, (q) = d 2 E;((\)/dq g . We notice that all such 
fluctuations arc correlated via the bandstructure of the device, since they originate 
from the same random variable x. 

We apply these ideas to the important case study of the Cooper pair box (CPB), 
a superconducting device described by the Hamiltonian 

Ho(q g ,J) = ^2(q g -n) 2 \(t>n)((j>n\- ^(l^n+t)(^i| + h.c.). (5) 


Here h := V/ n |D„) (6 n \ is the number of extra Cooper pairs in a metallic island. 
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The parameter J is the relative strength of the Josepson tunneling changing n —> 
nil. According to its value several different implementations of superconducting 
qubits from “charge qubits” to “transmons” [17J-2IJ have been demonstrated. The 
other parameter q g can be tuned by an external voltage. The CPB is operated by an ac 
gate voltage which is coupled to the charge, V = 2 eh, playing the role of the system 
operator. 

The Hamiltonian Ha(q rj , J) = JT E, i) (i\ is symmetric for charge parity trans¬ 
formations at half-integer and integer q g . Here the selection rule := (i\h \j) = 
[1 — (—1 y + i]riij holds, preventing pump coupling. On the other hand working 
at q g = 1/2 guarantees the maximal protection from charge noise because of the 
suppression of A,;(q) in Eq.([4]). Larger values of J 3 > 1 suppress asymmetries at 
q g 7 ^ 1 / 2 , ensuring protection in a larger region of the space of parameters, where 


however the pump coupling is suppressed [ 121 . 

The numerically calculated 1121 efficiency of STIRAP (Fig. [3ji) shows that a 
Lambda configuration allowing population transfer ~ 80% is achievable in a CPB 
with J = 1 by operating at a symmetry breaking bias q g / 1 /2, despite of the re¬ 
duced protection from low-frequency noise. In this regime, q g < 0.49, only linear 
fluctuation of detunings matter, i.e. linear terms in Eqs.Q are considered. We used 
HoT = 15 to guarantee good adiabaticity, with figures of noise and couplings consis¬ 
tent with measurements of the decoherence in the qubit of Ref. |22| . In this regime, 
Markovian emission |T4j 15 1 or absorption [12] channels are not effective, whereas 
spontaneous decay in the trapped subspace, characterized by T \, is tolerably small. 
Instead near q g = 1/2 pump coupling is small and it would require T 3 > T\, thereby 
decay suppresses the efficiency. 


4. STRATEGIES OF PROTECTION AGAINST NOISE 


We now analize the tradeoff between efficient coupling and decoherence when 
parity symmetry is broken. Our analysis leverages on the results of last section, 
that the main mechanism of efficiency loss (besides lack of adiabaticity) are low- 
frequency fluctuations of energy levels. As discussed in Secj3] these correspond to 
fluctuations of detunings and can be represented on the diagram of Fig. [2]by a curve. 
For vanishing nominal detunings it has parametric equations ([4]) and passes the origin 
for x = 0. Notice that in general in artificial atoms fluctuations of detunings produced 
by a specific noise source show specific correlations. In the linear noise regime 
5 p /6 = a, where a = ^(q )/A i (q), thereby the signs of the derivatives Aj(q) of 
the bandstructure determine whether detunings are correlated (e.g. critical current 


noise in flux qubits) or anticorrelated (e.g. charge noise in the CPB) [16]. It has been 


shown [TTJ that according to such correlations, three different kind of LZ patterns 
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Fig. 3 - (a) Efficiency of STIRAP in the CPB for J = 1, showing the effect of low-frequency noise 
alone and with additional Markovian noise. Field amplitudes, corresponding to Rabi oscillations f= 
2n • 600rad/s in the Quantronium qubit |22j operated at the symmetry point, yield a sufficiently large 
pump frequency Do — fffi n O2(9g)/ n oi(0.5) for values of q g which slightly break the symmetry. Flere 
we consider linear fluctuations of detunings (thick solid line), while adding their quadratic fluctuations 
(squares) and fluctuations of dtp in linear (dashed line) and quadratic approximation (dots) has no 
effect. Markovian noise determines a smaller further reduction of the efficiency (thin gray lines). We 
used a variance a x = 0.004 of p(x) and T\ = 500,1000ns (b) Figure of merit Eq.|6| for design and 
optimal symmetry breaking against low-frequency charge noise in the CPB. We consider the region 
where charge noise is linear since closer to the symmetry point q g = 1/2, STIRAP is suppressed by 
spontaneous decay. 


may occurr. They are shown in Fig. [2] (left panel) and labelled with (a) for 6 P /6>1, 
(b) for d p /6 < 0 and (c) for 0 < S p /5 < 1. The efficiency dependence on parameters 
turns out to be different in the three cases, this being the key to classify effects of 
low-frequency noise and the specific strategies to suppress them. 

Since both couplings and fluctuations depend on the solution of the eigenprob- 
lem of F/()(q), we may seek for an “optimal” set of values q such that the symmetry 
breaking yields enough pump coupling still keeping decoherence tolerable. Referring 
to Fig. |2jwe formulate this condition by defining a two-photon linewidth [231 <5i, as 
the interval of 5 where coherent transfer is appreciable, for any fixed combination of 
the other parameters (see Fig. [2]). For efficient STIRAP low-frequency noise must 

induce fluctuations of 5 with small enough variance, as ~ \JA/a/. + ! 2 Bfa/ < A, / 2 - 
The linewidth can be estimated by evaluating the impact of unwanted transi¬ 
tions between adiabatic states. In this way Vitanov et al. |23| found the scaling law 
5i ~ d(r )flo \JKp + Kg valid for 5 P = 0, and roughly holding on the lines 5 p /5 = a in 
the region (c) of Fig. |2j left. In the same way one can derive that <5i ~ d!(r, k)Qq k p 

in the region (b) of anticorrelated detunings [ 12], whereas <5i « d"(r, k s in the 

_ 2 

region (a) of correlated detunings. The dependence on k = k s / n p in the prefactors 
turns out to be weak. 
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In the case of CPB, since A\ and A 2 have different sign, charge noise de¬ 
termines anticorrelated fluctuations of detunings, and good transfer efficiencies are 
achieved for large values of the ratio 


5 Kpfio ^ 2 hq 2 (J,(Jg) 

as as \J + ( J ,q g )<4 


This is a figure of merit for STIRAP efficiency (see Fig. [ 3 J 5 ) which can be used for 
seeking optimization of both the design of the device and the symmetry breaking of 
the Hamiltonian modulated on-chip by the bias q g . 

The above analysis also suggests that effects of charge noise in a CPB can be 
minimized by increasing k p only. This is a specific way of decoupling dynamically 
noise sources, responsible for anticorrelated (5,d p ) fluctuations. Indeed it is clear 
from Fig. [2jo (the blue curve) that increasing n p the efficiency grows in the region (b), 
suppressing anticorrelated fluctuations of detunings. This happens because non ideal 
STIRAP occurs via LZ tunneling along the pattern (b) in Fig. |2} being suppressed 
for increasing 5 when the avoided crossing builts on, and being restored if the gap 
shrinks due to a larger il p . 

This analysis can be extended to the main different designs of superconduct¬ 
ing artificial atoms, and to each specific low-frequency noise source. These latter 
are classified according to the (6,6 P ) correlations they determine. For instance flux 
noise in flux qubits yields anticorrelated (5,6 P ), as for charge noise in the CPB, and 
increasing Q p yields dynamical decoupling. Instead critical current noise in CPB 
and flux qubit determine correlated (S, 5 P ) fluctuations, requiring larger Q s . In phase 
qubits both critical current and bias current noise yield correlated (6,6 P ) fluctuations 
dynamically suppressed by a large Q s . 

In real superconducting artificial atoms, where more than one noise source is 
present, the two strategies can be combined. Protection from noise producing an¬ 
ticorrelated ( 6 , 6 P ) fluctuations can be achieved by the optimal symmetry breaking 
strategy, since dynamical decoupling is limited by the insufficient coupling Q p . Pro¬ 
tection from noise producing correlated (5,5 P ) fluctuations can then be obtained in¬ 
creasing which is not limited by selection rules. 

It is easy to extend this analysis to artificial atoms driven in different field con¬ 
figurations. For instance for population transfer in the Ladder scheme(Fig.[l]) scheme 
one associates 5 (S p ) with the second (first) excited state, which allows to identify the 
relevant correlations between detunings. 
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5. IMPLICATIONS OF NON-MARKOVIANITY 

The picture of the last section relies on the non-Markovianity of BBCN. We 
remind that low-frequency noise is the main source of dephasing in artificial atoms. 
BBCN explains distinctive experimental observations in quantum bits | T6||22||24p5| . 
Moreover design of low-decoherence qubits relies on protection from non-Markovian 
noise. Both optimal tuning |26|27]| and dynamical decoupling [28 291 have been 
exploited for entangled states We generalize this ideas to protection of three-level 
coherence, obtaining a rich scenario. 

It is important to point out the different impact on STIRAP of non-Markovian 
dephasing, as discussed in this work, and Markovian pure dephasing as described 
by the standard Master Equation approach. This latter problem has been studied in 


Ref. [30], including only the dephasing rates 7 ^. For large enough OqT populations 
at the end of the protocol were found to be 


PnM = ^ e - 375lT2/8r 


Poo(00) = p 22 (00) = i - ie-W 2 / 8 " ( 7 ) 

i.e. dephasing determines efficiency losses which do not depend on the peak Rabi 
frequencies. Therefore Markovian dephasing cannot at all accout for the scenario 
presented in Sec. [4] 

In Fig. [4] we plot the final populations of the bare states comparing Markovian 
(pa) and non-Markovian (P t ) pure dephasing, in the entire relevant range of 0 (l . 
Noise produces in both cases the same qubit dephasing time T 2 , which is relatively 
large. For large Hq while for non-Markovian noise Pi is completely recovered, for 
Markovian noise it saturates to a smaller value given by Eq.(|7]). For small Qq the 
protocol fails in both cases, due to insufficient adiabaticity. 


6. CONCLUSIONS 

In this work we discussed effects of BBCN noise in three level artificial atoms. 
In particular we studied the tradeoff between protection from low-frequency noise, 
enforced by symmetries of the Hamiltonian, and the implied selection rules which 
are the main obstacle to the implementation of a Fambda scheme. Being based on 
two-photon absoiption and emission, the Fambda scheme allows performing tasks as 
transduction of photons in the pm/mm range. We have studied STIRAP since it is a 
benchmark advanced protocol. It is also the basis of other protocols from preparation 


of superpositions [ 15 1 to transfer of wavepackets and manipulation of photons, with 
still unexplored potentialities for quantum information and quantum control. 

We have shown that model for BBCN noise decoherence in the “trapped sub¬ 
space” span{|0), 11)} plays a major role, a conclusion which holds for all Fambda, 
Ladder and Vee schemes. Strategies to defeat noise in qubits can then be general- 
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Fig. 4 - Final populations of STIRAP with Markovian (pa) and non Markovian (Pj) noise. The for¬ 
mer is the solution of a Master Equation with rate 710 = I/T2, leading to exponential decay of qubit 
coherences. The non-Markovian noise is simulated taking a distribution of detunings corresponding to 
a x = leading to Gaussian decay with the same T 2 . 

ized to three-level systems. We presented two stetegies, namely optimal symmetry 
breaking and continuous dynamical decoupling, which can be integrated to minimize 
the effects of anticorrelated and correlated parametric fluctuations of the artificial 
atom bandstructure. Relying on non-Markovianity of BBCN, our results suggest that 
features of the scenario of STIRAP with BBCN, as the predictions on the peculiar 
dependence on control knobs described in Sec. [4j could be used to probe aspects of 
non-Markovianity of the solid-state evironment. 

Finally, we mention that artificial atoms allow for new unconventional schemes 
to manipulate a Lambda system, bypassing hardware constraints and allowing to per¬ 
form STIRAP at protected symmetry points and with always-on fields. The strategies 
to defeat noise presented here could be successfully applied also in these cases. 
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